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Abstract 

We present numerical solutions of the Gross-Pitaevskii equation corresponding to reconnecting 
vortex lines. We determine the separation of vortices as a function of time during the approach to 
reconnection, and study the formation of pyramidal vortex structures. Results are compared with 
analytical work and numerical studies based on the vortex filament method. 
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I. MOTIVATION 



The discrete nature of vorticity in quantum fluids (superfluid ^He, ^He and atomic Bose- 
Einstein condensates) makes it possible to give physical meaning to vortex reconnections. 
The importance of vortex reconnections in turbulence cannot be overestimated. Reconnec- 
tions randomise the velocity field thus sustaining the turbulent vortex tangle in a steady 
state [l| and changing the flow's topology jj]. By triggering a Kelvin wave cascade 3-7|, 
reconnections are responsible for the decay of turbulent kinetic energy at very low tempera- 
tures [8|. The dependence of the vortex reconnection frequency on the observed vortex line 

nn 

density is therefore an important quantity [9|, I10| . 



In quantum fluids, the existence of vortex reconnections was first 



1|; it was then theoretically demonstrated by Koplik and Levine 



conjectured by Schwarz 



observed at the University of Maryland 12|, |l3| using micron-size solid hydrogen tracers 



and experimentally 



trapped in superfluid vortex lines. More precisely, the Maryland investigators inferred the 
existence of reconnections by the sudden movement of trapped particles. In doing so they 
relied on numerical calculations by de Waele and Aarts [l^ , who reported that the distance 
between two vortices, 6, varies as 

S{t) = (K/27r)V2yt^, (1) 

where k is the quantum of circulation, t is time and to is the time of reconnection. This 
{tQ—ty^'^ scaling agrees with the approximate analytic solution of the nonlinear Schroedinger 
equation (also called the Gross-Pitaevskii equation, or GPE for short) found by Nazarenko 



and West [15|. De Waele and Aarts also claimed that, on their way toward a reconnection, 
the vortices form a universal pyramidal cusp structure which is independent of the initial 
condition. 



The numerical calculations of de Waele and Aarts [l4| were based on the vortex filament 
model pioneered by Schwarz which assumes that the vortex core is much smaller than 
any other length scale in the problem. The vortex filament model thus breaks down in the 
vicinity of a reconnection. It is therefore instructive to repeat the calculation seeking full 
numerical solutions of the GPE, which does not suffer such limitation and is a well-known 
model of superfluidity. This is indeed what we plan to do in this paper, together with 
revisiting the claim of universality of the reconnecting pyramidal cusp geometry. 
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II. MODEL 



The single-particle complex wavefunction ip(r,t) for bosons of mass m obeys the 3- 
dimensional time dependent GPE 

where r is the position, Vq is the strength of the delta function interaction between the 
bosons, and E is the single particle energy. The wavefunction is normalised by the condition 



that / \tp\'^dV = N. If we let r = h/{2E) be the unit of time, a = h/V2mE the unit of 



length, and ipo = y E/Vq the unit of ip^ we obtain the following dimensionless GPE 

-2^^ = VV + (1-|^I>. (3) 

We solve Eq. [3]in an (x, z) periodic box, using fourth order centred differences and Runge- 
Kutta-Fehlberg adaptive time stepping. -|-he typical numerical resolution is 100 x 300 x 120 
for example for the first calculation described in Sec. Illli More details of the numerical 



technique and the method used to set up initial conditions will be described elsewhere 16 |. 



III. RING RING RECONNECTIONS 



In a first set of numerical experiments, we consider the same planar vortex rings con- 

n 

figuration of de Waele and Aarts p^. Two vortex rings of (dimensionless) radius i? = 30 
are located on the z plane, initially separated by the (dimensionless) distance 5(0) = 20, 
and moving in the positive x direction. Fig. [1] shows that, during the evolution, the rings 
turn into each other. The pyramidal cusp which forms (and which eventually will lead to a 
vortex reconnection) trails behind the rest of the moving vortex configuration. 

By symmetry, the reconnection is on the 2; = plane, which facilitates the computation of 
the minimum distance between the vortex rings, (5(t), as a function of time t, until, at t = to, 
the rings reconnect. We repeat the calculation for different initial separations 5(0) = 30, 40, 
50 and 60. Fig. [2]shows a log-log plot of 5{t) vs t and confirms that indeed 5{t) ~ (to — tY^'^^ 
in agreement with Eq. [1] 

We fit our 5{t) vs t data to the form 
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FIG. 1: (Colour online). Ring-ring reconnection. Isosurfaces of = 0.7 (dimensionless units) 
showing two rings moving along the positive x direction at (dimensionless) times t = (left) and 
t = 81 right. Note the formation of the reconnecting cusp. 
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FIG. 2: (Colour online). Reconnections of vortex rings. Log-log plot of dimensionless distance 6{t) 
between the rings vs dimensionless time t, for various initial distances 5(0) = 20 (red crosses), 30 
(green diagonal crosses), 40 (blue asterisks), 50 (violet squares) and 60 (pale blue circles). 



5{t) = A^{to-t)[l + c{t^-t)l (4) 

where A and c are parameters. We find that A ^ c, and that A increases and c decreases 
for increasing 5(0). More precisely we obtain A = 2.66, 2.78, 3.18, 3.43 and 3.63, and 
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FIG. 3: (Colour online). Pyramidal reconnection cusp and definition of the angle (which a 
vortex makes with itself) and the angle 02 (which a vortex makes with the other vortex). 




FIG. 4: Reconnections of vortex rings. Reconnecting angles (pi (top, circles) and <j)2 (bottom, 
squares) vs initial distance (5(0). The lines are to guide the eye. 

c = 2 X 10"^ 1.06 X 10~^ 0.87 x 10^^ 0.67 x 10"^ and 0.52 x 10"^ respectively for 5(0) = 20, 
30, 40, 50 and 60. The errors are approximately ±0.08 for A and ±0.3 x 10"'^ for c, and 
arise mainly from the inaccuracy in localising the axes of the vortex lines. 



The Maryland group reported their results in the (dimensional) form 13 1 



Sditd) = AM\lt^d{tdQ - td)[l + CM{tdO - td)], (5) 

(where we introduced the subscript "d" to stress that a quantity is dimensional), where 
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FIG. 5: Initial vortex ring-line configuration to define the angle 6. 



A^i ^ 1.2 and cm ~ 0.5 are their fitting parameters. The numerical results of de Waele 
and Aarts jl4| correspond to Am = ~ 0.4. 

Since dimensional and dimensionless length, time and circulation are related hj 6 = S^/a, 
t = tfi/r and n = k^/ {(i'^/t) (where Kd = 10~^ cm^/s and n = 2tt are respectively the dimen- 
sional and dimensionless circulations), our (dimensionless) results agree with Maryland's 
(dimensional) results if A = AMyTK^/a'^ = AmV^tt ^ 3 and c = tcm ~ 10~^^ ~ 0, which 



is indeed the case (the value of E needed to get r to compare c with cm can be estimated 
from the definition of the healing length a and the fact that a ^ 10~^ cm). 

Following de Waele and Aarts [l^ , we define the intravortex angle 0i and the intervortex 
angle 02 of the pyramidal cusp formed by the two vortices as they approach the reconnection. 
These angles are computed from the dot products of the vector r2 with the vectors ri and 
r3, where ri = Pq — Pi, r2 = Pq — P2; and = Pq — P3 as shown in Fig. O If the total 
angle (ptot = 2(0i + 02) is equal to 360°, then the reconnection is "flat." If (ptot < 360° the 
vortices form a pyramid. 

Fig. m shows that reconnecting angles are approximately constant, 0i ~ 112° and 02 ~ 
61°, independently of the initial distance between the vortex rings, again in agreement with 
the finding of de Waele and Aarts [l^, although their numerical values are slightly different: 
their 0i increases from 115° to 135° over a change of two orders of magnitude of the initial 
distance, whereas 02 ~ 25°. 
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FIG. 6: (Colour online). Ring-line reconnections. Log-log plot of (dimensionless) distance 6{t) 
between the rings vs (dimensionless) time t, for initial angles 6 = 60° (red crosses), 75° (green 
diagonal crosses), 90° (blue asterixes), 105° (violet squares), 120° (light blue circles), 135° (yellow 
pointing up triangles), 150° (brown pointing down triangles) and 160° (orange diamonds). 

IV. RING LINE RECONNECTIONS 

In a second set of numerical experiments, we consider the interaction of a vortex line (set 
in the middle of the computational box and aligned in the z direction) with a vortex ring 
(located in the y,z plane, with centre at y = D) which moves in the positive x direction). 
The geometry of the initial configuration can be parametrised by the angle 6 between the 
vortex ring and the vortex line, see Fig. [Sj 6 ranges from 6 = 0° (the vortex ring passes to 
the left of the line) to ^ = 180° (the vortex ring passes to the right of the line). 

Fig. [6] shows that, for ring-line reconnections, the distance between the vortices is very 
similar to what is determined for ring-ring reconnections. Fitting as in Eq. H] we have 
A = 2.15, 2.28, 2.28, 2.02, 2.05, 2.01, 2.27 and 2.25, and c = -7.2 x 10-^ -7.2 x lO'^, 
-6.9 X 10~^ -7.2 X 10"3, -7.6 x 10"^ -8.9 x 10"^ -7.0 x 10"^ and -7.2 x 10"^ respectively 
for initial angles 9 = 60°, 75°, 90°, 105°, 120°, 135°, 150° and 160°. 

In contrast to what happens in ring-ring reconnections, there is less evidence for a uni- 
versal pyramidal cusp. Fig [7] shows that 0i and 02 depend on 6, that is to say, on the 
initial geometry. Fig [8] shows that the total angle (ptot (which would be 360° for a "flat" 
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FIG. 7: Ring-line reconnections. Reconnecting angles (pi (top, squares) and (p2 (bottom, circles) 
vs initial angle 9 (all angles in degrees). The lines are to guide the eye. 




FIG. 8: Ring-line reconnections. Total cusp angle (ptot vs initial angle (all angles in degrees). 
The line is to guide the eye. 

reconnection), is not constant. 

It is interesting to note that, if ^ < 50°, when the vortex ring is sufficiently close to the 
vortex line, the part of the vortex ring which is near the vortex and the part of the vortex 
line which is near the ring are almost parallel; since they have the same circulation, they 
tend to rotate about each other. This motion perturbs and deflects the vortex ring, which 
passes past the (perturbed) vortex line without any reconnection. Because of this curious 
effect, shown in Fig. [9], we have no data for 6 < 50° in Figs. [7] and [HI 



8 





(e)t=88 (f)t=110 



FIG. 9: (Colour online). Ring-line reconnection for 6 = 45°, seen from the bottom, looking up 
the z axis. Isosurfaces of IV'P = 0.7 (dimensionless units) at (domensionless) times t = 50, 70, 90 
and 110. The vortex ring (in the y,z plane) moves in the positive x direction past the vortex line 
(aligned along z) without reconnecting. 
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V. CONCLUSIONS 



By numerically finding nonlinear solutions of the governing GPE, we have confirmed 
that, as vortex filaments approach a reconnection, their distance scales as 5 ~ {to — t)^^"^, as 
predicted by Nazarenko and West 15| . By fitting our 6 vs t data to Eq. HI we find better 
quantitative agreement with Maryland's observations 13| than obtained with the vortex 
filament method fl^ . 

Finally, we find that reconnections involve the formation of a pyramidal cusp. However, 
on the length scales which we can explore using the GPE (which are different from length 
scales described by the vortex filament method), the angles of this cusp are not universal as 
previously claimed 1^ . 
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